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We consider a disordered two-dimensional system of independent lattice electrons in a perpendicular 
magnetic field with rigid confinement in one direction and generalized periodic boundary conditions 
(GPBC) in the other direction. The objects investigated numerically are the orbits in the plane 
spanned by the energy eigenvalues and the corresponding center of mass coordinate in the confined 
direction, parameterized by the phase characterizing the GPBC. The Kubo Hall conductivity is 
expressed in terms of the winding numbers of these orbits. For vanishing disorder the spectrum 
of the system consists of Harper bands with energy levels corresponding to the edge states within 
the band gaps. Disorder leads to broadening of the bands. For sufficiently large systems localized 
states occur in the band tails. We find that within the mobility gaps of bulk states the Diophantine 
equation determines the value of the Hall conductivity as known for systems with torus geometry 
(PBCs in both directions). Within the spectral bands of extended states the Hall conductivity 
fluctuates strongly. For sufficiently large systems the generic behavior of localization-delocalization 
transitions characteristic for the quantum Hall effect are recovered. 
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I. INTRODUCTION 

The two-dimensional electron gas in the presence of 
a strong perpendicular magnetic field provides a rich 
variety of physical phenomena. The most prominent 
is the quantum Hall effect MJ^. Another fascinating 
phenomenon is the self-similar energy spectrum (Hofs- 
tadter's butterfly) brought about by a periodic lattice 
potential and the magnetic field [^JJ] . These phenomena 
are due to two features originating from the magnetic 
field: chirality and a characteristic length scale. The chi- 
rality is due to the axial vector character of the magnetic 
field and is responsible for the Hall effect in the presence 
of an electric field and current carrying edge states in 
confined systems. The length scale is defined by the size 
of the cell associated with a single flux quantum and is 
called the magnetic length. In a periodic potential com- 
mensurability effects occur with respect to the ratio of 
lattice constant and magnetic length. In the presence 
of disorder also the localization length comes into play 
which measures the spatial extension of wave functions. 
The natural way to investigate the interplay of these dif- 
ferent length scales is to study the (dissipative and Hall) 
conductances of multi-terminal systems. Recently, a two- 
terminal conductance calculation for a clean system was 
performed by Skjanes et al. [|| and a multi-terminal con- 
ductance calculation has been performed by Aldea et al. 
|@|. For the time being the attempts of calculating di- 
rectly the multi-terminal conductances of disordered sys- 
tems are limited by computational practicability. There- 
fore, in the present paper we restrict ourselves to investi- 
gate the energy spectrum, its sensitivity to changing the 
boundary conditions and the Kubo Hall conductivity. 

The system to be studied is a two-dimensional An- 



derson tight-binding model with Pcicrls substitution (in 
Landau gauge) and generalized periodic boundary con- 
ditions (GPBC) in one direction and different bound- 
ary conditions in the other direction. This system has 
been investigated by many authors and several results 
are well established. In the limit of vanishing disorder the 
Schrodinger equation reduces to a one-dimensional finite 
difference equation known as Harper's equation which 
can be studied effectively by transfer matrix methods. 
Imposing periodic boundary conditions in both directions 
corresponds to a torus. If one of these is replaced by 
rigid wall confinement we are dealing with a cylinder. 
For both geometries the energy spectrum is periodic in 
the flux per unit lattice cell (denoted by a) with pe- 
riod 1. Concerning the torus the following results are 
established: (i) In the case of commensurability, i.e. ra- 
tional values a = p/q the energy spectrum consists of 
q Harper bands separated by gaps. The gap structure 
gives rise to Hofstadter's butterfly (ii) The energy 
spectrum is symmetric with respect to zero energy, (iii) 
The Hall conductivity is antisymmetric with respect to 
zero energy (referred to as particle- hole symmetry), (iv) 
Within the band gaps the Hall conductivity is quantized 
in integer multiples of e 2 /h\ its value is given by a topo- 
logical quantum number, the first Chern number on the 
magnetic Brillouin zone. The integer can be determined 
from the Diophantine equation Q 



pt n - qs n 



(1) 



Here t n , s n are integers corresponding to the n-th gap 
and the Hall conductivity there is an = t n (e 2 /h). The 
solution of Eq. (Q) is to be taken under the requirement 

|*n| < <Z/2. 

For the cylinder geometry the energy spectrum has 
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also been shown to consist of q Harper bands but the 
gaps are filled with states localized near the edges in the 
confined direction (x), and extended in the other. These 
edge states carry a chiral current and give rise to the 
same quantization of the Hall conductivity as obtained 
for the torus geometry || . Recently H an interpretation 
of the quantization in terms of topological invariants cor- 
responding to these edge states has been proposed and re- 
lated to the Chern numbers of the torus geometry. How- 
ever, this relation could be established for no disorder 
only (see also [[To[). 

So far, most studies have been focused on the ideal sys- 
tem where analytical results can be obtained (especially 
the Diophantinc equation for the Hall conductivity). It 
has been argued that disorder will not change the quan- 
tized value of the Hall conductivity as long as the disorder 
does not close the corresponding gap. This is certainly 
true for the torus geometry where the Hall conductivity 
is given by the topologically stable Chern number as long 
as the Fermi energy is within an energy gap. The argu- 
ment is however not proven for the cylinder geometry 
where no such gaps exist. Moreover, in the torus ge- 
ometry disorder leads to levels within the gaps, although 
the corresponding states are expected to be localized and 
so do not change the quantization as given by the Dio- 
phantine equation. For the torus geometry the following 
scenario is expected to occur: On increasing the disor- 
der strength localized states appear first in the tails of 
the disorder broadened Harper bands. When the disor- 
der is further increased the Harper bands begin to overlap 
strongly and form new bands with possibly less regions of 
extended states. Yet, for Fermi energies at which a small 
density of localized states is situated in a previous gap 
the topological stability of the Chern number still deter- 
mines the Hall conductivity via the Diophantinc equation 
(cf. Sec. VIII E in ||). 

For the disordered cylinder model (cf. jl2|,|l3|) we can- 
not rely on the topological stability of Chern numbers 
when discussing the Hall conductivity. Therefore, in the 
present work we take a different starting point: we inves- 
tigate the Hall conductivity in an appropriate version for 
the cylinder geometry (cf. Q ) and show that in this case 
also it has a topological meaning as it can be expressed 
in terms of a winding number. This winding number has 
an intuitive interpretation: it is the number of oriented 
windings of those orbits, in a diagram of energy eigenval- 
ues versus the center of mass coordinate, X, that connect 
opposite edge regions. The orbits are parameterized by 
a continuous quantum number which labels generalized 
periodic boundary conditions in the direction of current 
flow (y-direction). This winding number is similar to the 
Chern number when interpreted as the winding number 
of zeroes of the wavefunction with respect to the change 
of boundary conditions []l5 16 1. We study numerically 
the influence of increasing disorder on the validity of the 
Diophantine quantization. We also investigate the Hall 
conductivity within the Harper bands and show that the 
Hall conductivity strongly fluctuates as a function of the 



Fermi energy. For weak disorder and finite system sizes 
an unexpected collapse of the Hall conductivity occurs in 
the centers of the bands. This is, however, no surprise 
when consulting our version of the Kubo Hall conductiv- 
ity: As long as states in the band centers are extended 
and do not display large spatial fluctuations, X will be 
concentrated at the center of the system. Then the Hall 
conductivity stays close to zero. We expect that as soon 
as the disorder is strong enough in finite systems or ar- 
bitrary small but finite in infinite systems the collapse 
will disappear and (apart from strong fluctuations) the 
Hall conductivity effectively interpolates between subse- 
quent the quantized values. The reason for this expec- 
tation is that extended states (which in large systems 
are restricted to the energy band centers) show strong 
spatial fluctuations which are extremely sensitive to the 
variation of boundary conditions, i.e. the center coor- 
dinate X of such states can be shifted to an arbitrary 
position within the system by an appropriate change in 
the boundary conditions. To examine the spatial fluctua- 
tions of extended eigenstates and the universal properties 
of the localization-delocalization transition we calculate 
the scaling behavior of the distribution of wave function 
amplitudes of extended states. This (multifractal) anal- 
ysis (cf. Chap. 12 of Q) confirms our expectation and, 
furthermore, we find reasonable agreement with the scal- 
ing exponents obtained for other model systems of the 
integer quantum Hall effect. 

The paper is organized as fo llows: In Sec. |n] we intro- 
duce the model and in Sec. Ill the methods of calculating 
energy orbits and the Hall conductivity. The quantiza- 
tion of the Hall conductivity is determined by the topol- 
ogy of the energy orbits. In Sec. IV we present numerical 



results for the energy orbits corresponding to ideal sys- 
tems with cylinder geometry and, by comparing to the 
torus geometry, we demonstrate how edge states appear 
in the energy gaps of the bulk system. We then turn 
over to weak disorder and investigate its influence on the 
topology of energy orbits. The Hall conductivity is unaf- 
fected for Fermi energies situated in the bulk energy gaps 
and is still given by the Diophantine equation, Eq. (|l|). 
For Fermi energies situated in the bulk energy bands the 
Hall conductivity drops down and fluctuates around zero. 
To see whether this unexpected behavior is still valid for 
larger systems with strong disorder we present calcula- 
tions of energy orbits in such systems in Sec. [v| We 
discuss the fluctuation properties of the wave function 
for extended states in the weak and strong scattering 
regimes, respectively, and we conclude that for the latter 
the amplitude fluctuations are strong and display uni- 
versal characteristics of localization-delocalization tran- 
sitions in quantum Hall systems. The Hall conductiv- 
ity interpolates between adjacent plateau values as soon 
as the system size becomes much larger than the micro- 
scopic scales of the model but still exhibits strong meso- 
scopic fluctuations. Our conclusions are summarized in 
Sec.lvi. 
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II. THE MODEL 



Consider a two-dimensional (2D) square lattice with 



lattice constant a and extensions L x = M x a, L y 



M„a. 



With a set of orthonormal quantum states |m) localized 
at lattice sites m = (m x ,m y ) the Hamiltonian of the 
tight-binding model of independent electrons reads 



H = Y,£m\m) (m| + 



m) (m 



(2) 



(m,m') 



where (m, m'} denotes nearest neighbors. The site ener- 
gies £ m correspond to the (random) potential energy and 
the hopping matrix elements t m m / to the kinetic energy. 
The magnetic field is included in the kinetic energy by 
the Peierls substitution 



±27rzam a: 



if m y = m' y and m' x = m x ± 1 (3) 
if m x = n 



and m' = m y ± 1 . (4) 



Here a, restricted to < a < 1, is the number of 
flux quanta h/e per unit cell. The characteristic length 
associated to a is the magnetic length Ib defined by 
a = a 2 /2irl B . Requiring commensurability between the 
lattice constant and the magnetic length restricts a to 
rational values 



a 



p/q. 



(5) 



We impose generalized periodic boundary conditions 
on the eigenfunctions in y-direction 

M »i>(m x ,m y ) (6) 



ip(m x 



M y ) = e 



and Dirichlet boundary conditions in x-direction (cylin- 
der) 



V>(0, m y ) = ip(M x + 1, my) = 



(7) 



or (for comparison) periodic boundary conditions in re- 
direction (torus) 



tp(m x + M x ,m y ) = ip(m Xl m y ) . 



(8) 



The parameter -d e [0, 1) generalizes the periodicity of the 
wavefunction to be periodic "up to a phase". "&/M y can, 
however, also be interpreted as a Bloch quantum num- 
ber []l4f . Indeed, by transforming the eigenvalue problem 
Hip = Etf) according to the replacement 



iP(m x ,m y ) — > e 2 ^l M ^ m ^(m x ,m y ) 



(9) 



one arrives at an equivalent problem with strictly peri- 
odic boundary condition in y-direction. The transformed 



hopping matrix elements read for m x = m' x 

j- _ ±2-Ki[am w +S/My) 



= m' y ±l 
(10) 



A third way to interpret the parameter is by considering 
the system as the surface of a cylinder embedded in a 3D 
space. With y denoting the coordinate along the circum- 
ference of the cylinder, -d describes an Aharonov-Bohm 
flux (in units of the flux quantum) along the cylinder axis 
(cf. @). 



III. TOPOLOGY AND THE QUANTUM HALL 
EFFECT 



We adopt the view of considering a family of Hamil- 
tonians £T(i?) (characterized by Eqs. (^j,[|[l0]) ) and wave- 
functions that are periodic in y-direction. For each mem- 
ber of that family the velocity operator v y ($) can be ob- 
tained by taking the derivative of H (-#) with respect to 
the Bloch- momentum Tid/M yi 



My dH{ti) 



(ii) 



To model a disordered system we take independent 
random values for the site energies e m according to a 
symmetric box distribution on the interval [— V/2, V/2]. 

For vanishing disorder, V = 0, and zero magnetic field, 
a = 0, the Hamiltonian describes an electron moving in a 
2D crystal. The energy band is symmetric around E = 
and has the total width 8. On introducing the rational 
flux a = p/q and adopting toroidal boundary conditions 
(Eqs. (^|),(||) with $ — 0) commensurate with the flux 
lattice, i.e. M x = ql with some integer I, the spectrum 
consists of q magnetic subbands (Harper bands) sepa- 
rated by energy gaps (cf. || and Fig. 1 of Ref. [Q) and 
the statements (i) - (iv) listed in the Introduction apply. 

For V ^ the full 2D problem can be solved by diag- 
onalizing (numerically) the representing (sparse) Hamil- 
tonian matrix which yields the eigenvectors ipi(m; i9) and 
eigenvalues £/($) labeled by the discrete index I. In the 
present work we analyze the model by means of energy 
orbits to be introduced in the following. 

Defining the center of mass coordinate (in x-direction) 
of a given eigenstate by 



(12) 



allows us to calculate the Hall conductivity an = <j yx 
by using an appropriate version of the Kubo formula 
(cf. Eq.(4.76)) in M) which, for zero temperature, reads 



d<d 



(13) 



where £p is the Fermi energy. This expression tells that 
once we know the "energy-orbits" (A;($), £i($)) we can 
read off the Hall conductivity at a given Fermi energy by 
determining the intersection points, (IW.eW = ep), of 
the orbits with the line of constant Fermi energy. 



oh 



E 



sign(i) X^ 



M~ 



(14) 



This formula tells that one has to add up the center of 
mass coordinates of the intersection points with a signum 
± depending on crossing from below (— ) or above (+) 
and finally divides by the width M x resulting in the 



3 



Hall conductivity (in atomic units). Note that the or- 
bits are (generically) closed, since level crossings occur 
with zero probabilit y in a disordered system (Wigner-von 
Neumann theorem ||L7| ). 

Those orbits which connect both edge regions of the 
cylinder are of particular importance: they contribute 
(almost) integer numbers to the Hall conductivity. 

These integer numbers are given by the winding number 
of orbits connecting both edge regions of the cylinder. The 
chirality of the orbits defines the sign of the integer, (+) 
for clockwise and (—) for counterclockwise windings. 

The slight difference between the center of mass co- 
ordinate of an edge state X and the geometrical edge 
(m x = 0, M x ) gives rise to corrections. However, they 
are of order (Ib/L x ) and furthermore they turn out to be 
artificially caused by using conductivity instead of con- 
ductance (cf. the discussion in E3] and Sec. 4.8 of p|). 
In our numerical calculations these corrections are still 



visible for system sizes of order 100a. Equation (13) has 
been analyzed previously in the context of continuous 
electron systems on the surface of a cylinder p4] , ^9| , ^8 |. 
We mention that the very fact of level anti-crossing leads 
to the presence of energy gaps (mini gaps) in the regime 
of edge states. Within these mini gaps the Hall conduc- 
tivity vanishes. A similar observation in a two terminal 
conductance calculation was denoted in [^) as "cracks in 
the wings' of the Hofstadter butterfly. However, these 
gaps are exponentially small in the ratio of magnetic 
length vs. system size (cf. |l^]) and will be ignored in 
the following discussions. 

We like to comment further on the topological quanti- 
zation of the Hall conductivity as expressed by the wind- 
ing numbers of energy orbits. A condition for the occur- 
rence of plateaus in the Hall conductivity is the existence 
of localized states in the bulk of the system. Such states 
are characterized by an exponential decay in space far 
off their center position. The decay length is called the 
localization length £. Localized bulk states correspond 
to orbits degenerated to points since the sensitivity of 
the center of mass to a change in boundary conditions 
becomes exponentially weak in the ratio £,/L x . If, for a 
given Fermi energy, all bulk states of the system are local- 
ized but extended edge states do exist, the winding num- 
ber of these edge states determine the quantized value 
of the Hall conductivity. This discussion is somewhat 
simplified since localized bulk states and extended edge 
states are not degenerate in energy. Rather, the states 
at a given energy are hybrids of both types of states (see 
Fig. ||a). As already mentioned in the Introduction the 
topological quantization in terms of winding numbers is 
in close analogy to the topological quantization in terms 
of Chern numbers occurring in the torus geometry. The 
Chern numbers are given by the winding of the zeroes 
of wave functions around the torus when two parame- 
ters (Bloch quantum numbers) for generalized periodic 
boundary conditions are changed [ p5|p^ ]: Moving the 
Bloch quantum numbers over their entire Brillouin zone 



(being also a torus) leads to an integer covering num- 
ber on the torus in real space. This integer is just the 
Chern number determining the Hall conductivity. Thus, 
the Hall conductivity counts the number of times the ze- 
roes of a wave function cover the torus when changing 
the boundary conditions. In the cylinder geometry the 
Hall conductivity counts the number of times the center 
of mass of a wave function covers the cylinder when the 
boundary conditions are changed. This analogy demon- 
strates the universality of topological quantization argu- 
ments for the Hall conductivity; they are not restricted 
to a particular geometry although for different geome- 
tries the topological quantization can manifest itself in 
different objects. 

A further indication for this universality comes from 
the observation that the distinction between localized 
and extended states can be based on topological prop- 
erties of either model. In the torus geometry localized 
states are characterized by a vanishing Chern number 
corresponding to a vanishing covering number for the ze- 
roes of the wave function. In the cylinder geometry lo- 
calized states are characterized by a vanishing covering 
number for the center of mass of the wave function. In 
a quasi-classical treatment (see e.g. |?0j) localized states 
are characterized by closed equipotcntial lines that do 
not wind around the cylinder. 

Equation ( |l3| ) not only helps in demonstrating the 
quantization of the Hall conductivity in a regime of lo- 
calization but is also suitable for discussing the Hall 
conductivity in the transition regimes between adjacent 
plateaus. Here complicated orbits lead to a strongly fluc- 
tuating Hall conductivity, as will be discussed in the next 
sections. It is worth mentioning that similar fluctuations 
are expected to occur in the torus geometry whenever the 
Fermi energy crosses a band of extended states (cf . ) . 

A quantity that is simpler than the Hall conductivity 
and usually still allows to distinguish between localized 
and extended states is the so-called Thouless number G 
(cf. (2^]) defined as the ratio of a typical energy shift 5e 
within a band and the average level spacing A 



G 



6s 
A 



(15) 



This quantity is similar to a dissipative conductance (in 
atomic units) since it becomes small compared to 1 for 
localized states. However, as we will see, for the cylin- 
der model there are extended edge states even in energy 
ranges where the bulk states are localized. Thus, the 
Thouless number will show some non-zero peaks indicat- 
ing that it is not a suitable criterion for a bulk state to 
be extended. Therefore we will focus in the following 
sections on the Hall conductivity, Eq. (|l3|). 



IV. WEAK DISORDER 

We begin our discussion with the limit of vanishing dis- 
order and focus on a significant difference between the 
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a) 



FIG. 1. Energy spectrum as a function of flux. 




FIG. 2. Energy branches as a function of wave number for 
flux a = 1/5. 



torus and the cylinder geometry. As has been pointed 
out in H the difference with respect to the spectrum is 
that in the cylinder model energy eigenvalues occur in 
the Harper gaps of the torus model. This can be seen in 
Fig. |l| where the band structure for the cylinder model is 
shown as a function of a (the number of flux quanta per 
unit cell). Regions with dense set of points correspond 
to the Harper bands while those with lines of points cor- 
respond to the gaps (see also Fig. 1 in Q and Fig. 4 in 
||). For vanishing disorder the translational invariance 
in y-direction yields a wave number k labeling the energy 
eigenvalues (cf. S). The existence of energy branches 
within the Harpergaps can be easily observed by plot- 
ting Si(k). In Fig. gsuch plot is shown for a — 1/5, d = 
and M x = M y — 17. One observes both the 5 Harper 
bands and some energy branches within the Harper gaps. 
To see that the existence of energy eigenvalues in the 




b) 



FIG. 3. Energy orbits for a system with no disorder and 
flux a = 1/5 for a torus geometry (a) and for a cylinder 
geometry (b). 



Harper gaps is due to edge states it is instructive to plot 
the energy eigenvalues versus the corresponding center of 
mass coordinate (in ir-direction) for a number of $ val- 
ues. The corresponding plot for the torus geometry is 
shown in part a) of Fig. |2|. This plot not only shows the 
Harper bands but also demonstrates that the center of 
mass, X, is distributed within the bulk of the system. 
For a system with no disorder the wavefunctions are pe- 
riodic in x-dircction j23|. Thus, X is distributed within a 
range smaller than the period. On introducing the cylin- 
der boundary conditions the picture changes. As can be 
seen in part b) of Fig. [|, the states which have energies 
within the Harper gaps are concentrated along the edges 
(see also Fig. 6 of Ref. ||). Changing d from to 1 can 
move an edge state from one edge to the opposite edge 
and, finally, back to the starting point. On the other 
hand, the states with energies within the Harper bands 
are less effected by the presence of an edge. 

The existence of edge states can be qualitatively under- 
stood by considering the Dirichlet boundary conditions 
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FIG. 4. Energy orbits for a cylindrical system with no dis- 
order and flux a — 3/7. 



hi 

b 




FIG. 5. Hall conductivity as a function of energy for disor- 
der strength V = 0.9, flux a = 1/9 and system size M = 20. 



together with the chirality inherent in the phase of the 
hopping elements. For an electron that starts to propa- 
gate close to an edge these items lead to a directedness 
of the most probable path along the edges. 

The energy orbits (£;($), -X"z($)) provide a very helpful 
picture of the spectral properties of the cylinder model 
and allow to calculate the Hall conductivity rather easily, 
as discussed after Eq. ( |l3| ) . For example the Hall conduc- 
tivity for a = 1/5 is 1, 2, —2, —1 in the first, second, third 
and fourth gap (from below) , respectively. These findings 
are consistent with the Diophantine equation ([!]). For 
values a = 1/q, with q large, the low lying Harper bands 
strongly resemble Landau bands of continuous electron 
models [ p4[ and the quantization of ajj in increasing inte- 
ger steps is no surprise. However, in the lattice model we 
can take a rational value with a dominator different from 
1, e.g. a = 3/7. Here the Diophantine equation predicts 
a sequence of an as follows: uu = —2, 3, 1, —1, —3, 2. As 
can be seen in Fig. ^ the result is still valid in the cylinder 
model and coincides with counting the winding numbers 
of orbits connecting edge states. In fact, in all the cases 
we analyzed the chirality of the orbits leads to the sign as 
predicted by the Diophantine equation. Thus, we recover 
the result obtained in |J that the Hall conductivity for 
the cylinder model with vanishing disorder coincides with 
the Hall conductivity for the torus model with vanishing 
disorder, as expressed by Eq. (|l]). The point we wanted 
to make here was that the concept of winding numbers of 
energy orbits is suitable to determine the Hall conductiv- 
ity. We will use this concept further on when increasing 
successively the strength of disorder. 

Let us first consider the case of weak disorder. By 
this we mean that the disorder strength V is small as 
compared to the strength of the hopping matrix element 
( set 1). Since for any given matrix size N the num- 
ber of diagonal elements and the number of non-zero 
off-diagonal hopping elements are both of order N, sig- 



nificant changes in global spectral properties (such as 
the density of states) occur when V becomes of order 
1. Note, however, that this definition of weak disorder 
does not apply to disorder related length scales. For ex- 
ample, once the system is large enough localized states 
will show up in regimes of low density of states. Thus, 
arbitrarily small amount of disorder can lead to finite lo- 
calization lengths £. From this point of view disorder is 
always strong. 

In the following we will concentrate on a cylinder sys- 
tem with a = 1/9 and M x = M y = 20. For V = 0.9 the 
Harper bands become broader but still do not overlap. 
In Fig. H the Hall conductivity is shown as a function of 
the Fermi energy for V = 0.9. The quantization of an 
in the gaps remains valid. However, the behavior of <jh 
within the broadened Harper band is striking: it drops 
down and fluctuates around zero. A similar behavior can 
be concluded from the plot of energy orbits for V = 0.6 
displayed in Fig. ||. Again the quantization in the region 
of bulk states energy gaps according to the Diophantine 
equation remains valid (an — e 2 //i(l, 2, 3, — 3, — 2, — 3)) 
but the variation of center coordinates X for extended 
bulk states is restricted to a range of approximately 10 
lattice spacings around the center of the system. There- 
fore, the Hall conductivity will not exceed the value 0.5 
within the band. The question arises if this behavior 
will continue for larger systems. If it does, the Hall con- 
ductivity will not interpolate between adjacent plateaus 
but drop down in the transition regimes. To answer this 
question we first have to understand the origin of the 
drop in the Hall conductivity. Equation ( |T3| ) relates the 
Hall conductivity to the range of center coordinates X 
for a given energy orbit. For systems of no disorder the 
amplitude of a wave function is periodic over the flux 
lattice. Consequently, the range of X is restricted to be 
less than the period p3j. As soon as we introduce dis- 
order we have to ask about the nature of extended bulk 
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FIG. 6. Energy orbits for disorder strength V = 0.6, flux 
a = 1/9 and system size M = 20 



states. It is clear from Eq. (13) that the range of X values 



can only increase provided the amplitudes of an extended 
wave function are strongly inhomogeneous. Only then an 
appropriate shift in the boundary condition parameter $ 
can lead to a drastic shift in the position X of the center 
of mass. 

The localization length £ of a wavefunction is the rel- 
evant length scale introduced by disorder. It is expected 
that for sufficiently large system sizes extended states 
can only exist in the center of the broadened Harper 
bands. All other bulk states are localized and do not 
contribute to er#. The bulk extended states are charac- 
terized by their localization length £ being larger than 
the actual system size. In the thermodynamic limit 
M = M x = My — > co the spectral width of extended bulk 
states is expected to shrink to zero as M~ x l v where v is 
the critical exponent of the localization length (cf. 
As long as the system size is so small that we hardly find 
localized states, even in band tails, the system size has to 
be considered as "microscopic" . The plot of the squared 
amplitudes (density plot) of two wave functions corre- 
sponding to the tail of the first band and to the center 
of the second band of bulk states are shown in Fig. 0a) 
and Fig. Qb), respectively. The first state shows a ten- 
dency to localize, however the localization length is of 
the order of the system size. The second state is delocal- 
ized and already displays some inhomogeneous amplitude 
fluctuations. However, for the system size of M — 20 the 
center of mass cannot be shifted over the entire scale 
of ^-coordinates. Nevertheless, the second state already 
indicates that the extended bulk states can become spa- 
tially inhomogeneous. 



b) 

FIG. 7. Density plot of two wave functions with energy in 
the tail of the first Harper band (a) and within the center of 
the second Harper band for the same system parameters as 
in Fig. 6. 



V. STRONG DISORDER 

To see the influence of stronger disorder on the 
Hall conductivity one can either consider larger systems 
with fixed disorder strength V or increase the disorder 
strength V. In both cases a reasonable amount of local- 
ized states show up in the low density of states regime. 
We changed the disorder strength and system size for 
a = 1/9 to V = 2 and M = 100, respectively The 
influence of the increased disorder strength can be ob- 
served in the density of states which is shown in Fig. ||. 
The Harper bands strongly overlap and the resolution of 
only the first and the last Harper gap is possible. These 
gaps are expected to be filled with localized states and 
form mobility gaps rather than spectral gaps. We will 
concentrate on the wave function properties for energies 
situated within the lowest mobility gap and in the center 
of the second Harper band. 

We find that within the gap all bulk states are now 
localized and only extended edge states appear, which 
meets the expectations formulated above. This can be 
seen in Fig. ||a) where the density plot of a wave func- 
tion corresponding to a gap energy is displayed. The plot 
reflects that the wave function has a bulk contribution 
concentrated on an area with diameter of approximately 
16 lattice spacings and an edge contribution which is ex- 
tended along one edge in y-direction. In contrast to the 
edge contribution the bulk contribution is almost insensi- 
tive to a change in the boundary conditions. For a certain 
value of the parameter -d the edge contribution will sweep 
to the opposite edge, resulting in a quantized Hall effect 
for Fermi energies situated in the mobility gap. The cor- 
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FIG. 8. Density of states for flux a = 
strength V — 2 and system size M = 100. 



1/9 and disorder 



responding energy orbits are shown in Fig. ^b) where the 
movement of edge states can be seen. 

The density plot of Fig. |lO|a) shows a state within the 
center of the second Harper band. It is a typical example 
of a critical eigenstate which shows strong spatial ampli- 
tude fluctuations. As a consequence the position of its 
center of mass X can be shifted over a large range of val- 
ues by an appropriate shift in the parameter as seen 
in Fig. 10b). By this observation we are led to conjec- 




FIG. 9. Density plot of a wave function with energy in a 
mobility gap (a) and the corresponding energy orbits (b). The 
system parameters are the same as in Fig. 8. 



b) 
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FIG. 10. Density plot of a wave function with energy in the 
center of the second Harper band (a) and the corresponding 
energy orbits (b). The system parameters are the same as in 
Fig. 9. 



ture that the corresponding Hall conductivity will inter- 
polate between adjacent plateau values when increasing 
further the system size. Still, the Hall conductivity shows 
strong (mesoscopic) fluctuations of 0(e 2 /h) related to 
the fluctuations of the center of mass coordinate X. For 
this conclusion it is essential that extended bulk states 
show strong spatial fluctuations. Such fluctuations are 
a generic feature of electron states at the localization- 
delocalization transition. In addition, these fluctuations 
appear on all length scales between microscopic scales 
and the (diverging) macroscopic localization length and 
consequently the corresponding density plot represents a 
multifractal (cf. Chap. 12 of ||). An important feature 
of critical states is that the typical density, defined as a 
geometric mean, P typ = exp (in |^| 2 ) scales in a universal 
manner with the system size M 



P tyP oc M~ 



(16) 



where ao is a critical exponent characteristic of the 
localization-delocalization transition at hand. The fact 
that ao is larger than the space dimension indicates the 
difference between average and typical values of densities 
and is a criterion for multifractality. 

In the work by Huckestein et al. |25|] the critical eigen- 
states of the lowest Harper band for a = 1/8 have been 
analyzed for a torus geometry. They found the multifrac- 
tal exponent ao to be consistent with values reported 
for a number of models related to the quantum Hall 
effect [^6|. A recent calculation ]27j yields a value of 
ao = 2.28 ± 0.02. For a few of the critical states in the 



two lowest Harper bands (a = 1/9, V = 2, M = 100) 
we calculated the exponent ocq and found values of ao 
compatible with those reported previously. This indi- 
cates that in the center of the Harper bands localization- 
delocalization transitions occur in accordance with the 
quantum Hall universality class. 



VI. CONCLUSION 

We presented numerical calculations of energy eigen- 
values and wave functions for disordered lattice electron 
systems on a cylinder surface in the presence of ratio- 
nal magnetic flux. We analyzed the Hall conductivity by 
means of closed orbits in the energy vs. center of mass 
plane. The orbits are parameterized by an Aharonov- 
Bohm flux along the cylinder axis. The spectrum is 
found to consist of disorder broadened Harper bands with 
Harper gaps filled by localized bulk states and/or delo- 
calized edge states. In the regime of these bulk mobil- 
ity gaps the Hall conductivity is quantized and given by 
the Diophantine equation, Eq. (|l|), originally derived for 
clean systems with torus geometry. 

In close analogy to systems with torus geometry the 
quantized Hall conductivity for systems with cylinder ge- 
ometry is shown to be related to topological quantum 
numbers. Whereas for the torus these are Chern num- 
bers, they are the winding numbers of the above defined 
orbits in the case of the cylinder. Thus, concerning the of- 
ten discussed "topological explanation" of the quantum 
Hall effect both geometries are shown to be essentially 
equivalent. 

Within the broadened Harper bands a regime of ex- 
tended bulk states exist. Regarding the drop in the Hall 
conductivity for very small and weakly disordered sys- 
tems, we do not expect that this drop will occur in a real 
multi-terminal measurement. Let us stress that the Hall 
conductivity assumes a uniform electric field within the 
system while a true conductance must be based on the 
self-consistent field within the system. As discussed in 
Ref. [U the difference between conductance and conduc- 
tivity is negligeble (0(Ib/L)) for discussing the plateau 
region, however this is not the case in the region of bulk- 
extended states. These bulk states are, for sufficiently 
large enough systems, multifractal with universal charac- 
teristics of critical states in quantum Hall systems. The 
multifractality is directly connected to the mesoscopic 
fluctuations of the Hall conductivity within the regime 
of critical states; the sensitivity of the center of mass to 
the Aharonov-Bohm flux leads to drastic fluctuations of 
the Hall conductivity. These fluctuations are of order 
e 2 /h and the distribution of the Hall conductivity in the 
transition regime between Hall plateaus deserves further 
studies. Although the very fact of mesoscopic fluctua- 
tions in the conductivity points to a similar effect for the 
conductance, a direct quantitative comparison is not pos- 
sible. The work by Aldea et al. , where multi-terminal 



conductances were calculated (for systems of rs 600 lat- 
tice sites), gives support to the expectation that a true 
multi-terminal Hall conductance fluctuates between sub- 
sequent plateaus. The experimentally observed flat dis- 
tribution pa ] of two-terminal conductances in the transi- 
tion regime is compatible with the fluctuations we found 
in our calculations of the Hall conductivity. They have 
been reproduced for the two-terminal conductance within 
a (one Landau band) network model of a quantum Hall 
system p9| . In our model similar fluctuations correspond 
to the strong amplitude fluctuations in critical (multifrac- 
tal) eigenstates. 
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